Thermal evolution of the full three-dimensional magnetic excitations in the multiferroic BiFe0 3 
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The idea of embedding and transmitting information 
within the fluctuations of the magnetic moments of spins 
(spin waves) has been recently proposed^ and experimen- 
tally tested* 2 - The coherence of spin waves, which describes 
how well defined these excitations are, is of course vital to 
this process, and the most significant factor that affects the 
spin-wave coherence is temperature. Here we present neu- 
tron inelastic scattering measurements of the full three- 
dimensional spin-wave dispersion in BiFeOa, which is one 
of the most promising functional multiferroic material 3 , 
for temperatures from 5 K to 700 K. Despite the presence 
of strong electromagnetic coupling, the magnetic excita- 
tions behave like conventional magnons over all parts of 
the Brillouin zone. At low temperature the spin-waves 
are well-defined coherent modes, described by a classical 
model for a G-type antiferromagnet. A spin-wave veloc- 
ity softening is already present at room temperature, and 
more pronounced damping occurs as the magnetic order- 
ing temperature T N ~ 640 K is approached. In addition, 
a strong hybridization of the Fe 3d and O 2p states is found 
to modify the distribution of the spin-wave spectral weight 
significantly, which implies that the spins are not restricted 
to the Fe atomic sites as previously believed. 

Coupling between static magnetic and ferroelectric orders 
in multiferroic materials has attracted tremendous recent in- 
terest because of its potential use in device applications^— If 
such coupling were to extend to dynamical properties, e. g. be- 
tween lattice vibrations and spin fluctuations, then a new type 
of excitation may emerge. Such hybrid excitations, or "elec- 
tromagnons," have already been studied in a number of multi- 
ferroic systems using dielectric and Raman measurements^— 
An even more fascinating aspect of this novel concept is the 
potential ability to control the spin waves using the electric 
fields generated by low-power dissipating electronic circuits 
instead of magnetic fields, a possibility that has been sug- 
gested for multiferroic BiFeOa at room temperaturer^*!! mak- 
ing this material of particular interest 3 -. Unfortunately optical 
measurements such as Raman scattering are only sensitive to 
long-wavelength (q ~ 0) excitations. Neutron scattering tech- 
niques, on the other hand, provides a powerful tool to probe 
these hybrid excitations over a large momentum space. There 
have already been some pioneering neutron scattering mea- 
surements on the magnetic excitations from powder BiFeC>3 
samples™, as well as single crsytal sample s 13 ' 14 but mea- 
sured only around room temperature. A more complete pic- 



ture of these modified spin-wave excitations, i. e. the electro- 
magnons, at all wave vectors Q, and especially their evolution 
with temperature is highly desired. 

BiFe03 becomes ferroelectric at Tc ~ HOOK. Below Tc 
this material exhibits an i?3 C 15 ' 16 crystal structure, which can 
be viewed as a (weakly) rhombohedrally-distorted perovskite 
structure in which the oxygen octahedra are also distorted. 
For the purpose of our magnetic analysis this is well approxi- 
mated by the pseudo-cubic unit cell shown in Fig.Q](a). The 
Neel temperature for BiFeOs is Tjy ~ 640 K, and the ordered 
magnetic phase has a cycloid modulation- with a long peri- 
odicity that is superimposed on a simple G-type antiferromag- 
netic structure. The magnetic Bragg peaks appear close to the 
magnetic zone-centers of Qaf = (0.5, 0.5, 0.5) (r.l.u.). 

In Fig. |2] (a)-(d), the magnetic excitations measured at 
T=5K along [001], [110], [111], and [210] directions are 
plotted. At base temperature, the spin-wave excitations are 
well defined, and disperse quickly when moving away from 
Qaf- One will see that although the features at higher ener- 
gies near the top of the bands are very well resolved, fine fea- 
tures at low energies for Q near Qaf cannot be accurately de- 
termined, due to our relatively coarse energy and wave-vector 
resolution. We then choose to model the magnetic Hamilto- 
nian by ignoring the effect of the long-period cycloid modu- 
lation, or any single-ion spin anisotropy in this system, which 
would only affec t 13,14 the magnetic excitations at low ener- 
gies and in a narrow range around Qaf ± 5, where 5 ~ 0.004 
along (110) are the modulation wave-vectors for the cycloid 
structure 18 . The Heisenberg spin Hamiltonian of the system 
therefore can be simplified to that of a classical G-type AFM, 
H = J J2 n m ^ n ' where we only consider the nearest 
neighbor interaction Ji, second nearest neighbor interaction 
J2, and the third nearest neighbor interaction J3, as shown 
in Fig. Q] (a). The interaction that dominates is clearly the 
nearest-neighbor exchange interaction Ji between Fe spins 
separated along the (100) direction, whose magnitude can 
be roughly estimated from the top of the magnetic dispersion 
band, while J2 and J3 can also be obtained if the full disper- 
sion is known (see Methods section for details). We have per- 
formed an overall fitting to our data using the model described 
in the Methods section to determine the exchange constants. 
The model calculations are shown in Fig. |2](e)-(h), with fit- 
ting parameters shown in Table U The calculated dispersion 
curves are also plotted in Fig. [2] as solid lines. It is apparent 
that our model provides a good description to the data at 5 K. 
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FIG. 1. (Color online) Schematics of the BiFeC>3 crystal structure (a), and spin densities from conventional Fe ionic model (b) and from first 
principles Wannier analysis taking into account the O 2p orbital (c). 



The fact that Jj_, J 2 and J 3 are all positive, suggests that there 
is frustration in the system. Specifically, the face-diagonal Fe- 
Fe spins are aligned along the same direction in spite of the 
fact that the interactions between them are antiferromagnetic. 

Upon heating the spin-wave is modified gradually. The in- 
tensity of low energy magnetic excitations is boosted by the 
Bose factor from 5 K to 140 K and then 300 K, while the over- 
all shape still remain unchanged [see Fig. [3(a) - (c)]. This is 
expected and can be understood as intensities from the mag- 
netic Bragg peak are shifted into the inelastic channels when 
the static order becomes less robust with heating. At higher 
energies near the band top, the change is minimal between 
5 K and 140 K, but at 300 K the band top is noticeably lower 
[Fig. |3](c)], and the excitations become broader in energy. The 
softening of the spin-wave modes is due to the renormaliza- 
tion that occurs at non-zero temperatures, which is essentially 
a result of interactions between the excited magnons, and can 
be modeled with various theories. Instead of going into details 
of the renormalization, we would still use the same simple 
spin-wave model that describes the low-T data well, but with 
a set of renormalized effective exchange parameters shown in 
Table U At 140 K the effective exchange parameters are virtu- 
ally the same as those obtained at 5 K. The effective exchange 
constants J\ and J2 both show a clear drop at T=300 K, and 
a broadening in energy of Y ~ 5 meV is required to obtain a 
good fit to the data. Our results indicate that although the spin- 
waves are still relatively well defined at room temperature, the 
lifetime (inverse of the energy width) of the spin-wave mode 
does decrease significantly compared to that at base tempera- 



ture. The magnitude of the spin- wave softening is smaller than 
the softening observed of the electromagnon modes 8 by Ra- 
man measurements near q = 0, but note that the energy of the 
q = mode in the Raman work is also related to the cycloid 
wave-vector— 6, which also drops with heating-^. The clear 
spin-wave velocity softening at 300 K makes the low tempera- 
ture measurements important to obtain the ground state values 
of the exchange constants and related parameters. 

With further heating towards TV ~ 640 K, the spin-waves 
become heavily damped, as shown in Fig.|3](d) and (e). The 
data taken at 580 K already show that although the scattered 
intensities are more intense around Qaf due to the dominat- 
ing AFM correlations, well defined spin-waves are no longer 
present. While the Neel temperature marks the disappearance 
of static magnetic order in the system, the coherence of spin- 
wave excitations already diminishes at temperatures well be- 
low Tjv. 

In addition to probing the coherence of the spin-wave ex- 
citations, which are collective magnetic responses from the 
system, our neutron scattering measurements also provide in- 
formation on the local spin structure around each Fe atom. 
The magnetic form factor |F(Q)| 2 describes the Fourier trans- 
formation of the spatial distribution of spin (density) around 
each magnetic ion (in our case, Fe 3+ ), into momentum space. 
In Fig.|4](b-d), we show the magnetic form factor for a Fe 3+ 
ion with the blue solid lines. This magnetic form factor is 
isotropic and only depends on the length of Q. However, re- 
cent analysis of inelastic neutron scattering data from high-Tc 
cuprates suggests^ that the spin distribution around magnetic 
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FIG. 2. (Color online) Magnetic excitations from BiFeC>3 at T=5 K. 
The vertical axes denote energy transfer, while the horizontal axes 
are Q - Qaf for Q along [100], [1 10], [1 1 1], and [210] directions 
for (a), (b), (c) and (d), respectively. The strong background appear- 
ing at large Q comes from phonon scattering, (e) to (h) are model 
calculations described in the text. 



ions in metal oxides can be greatly affected by the presence 
of the oxygen ion, and could result in different effective mag- 
netic form factors. 

To compute the magnetic form-factor that includes the ef- 
fects of the crystalline environment, a first principles Wan- 
nier analysi o 20 ' 21 is employed. From LDA+U calculations of 
BiFeO ^ 22 ' 23 the Fe is found to be in the high-spin 5 = 5/2 
state. The magnetic form-factor is obtained by Fourier trans- 
forming the spin density of the low-energy Wannier functions 
of Fe. In Fig. Q] (b) and (c), we compare the atomic spin- 
density with the Wannier function derived spin-density. The 
main difference is that the Wannier function derived spin- 
density spreads into the neighboring oxygen atoms, which 
suppresses the magnetic form-factor at none-zero momenta. 
Direct comparison of the form factors are plotted in Fig.|4](b)- 
(d). We can see that for most Q, the Wannier form factor is 
smaller than the previously used ionic form factor. An exper- 
imental verification of the new form factor is shown in Fig. [4] 
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where 



(a). Here we plot the ratio R(Q) = I{Qo)/lF (Qo) | 
/ (Q) is magnetic scattering intensity measured at wave-vector 
Q, while Qq is the wave-vector equivalent to Q in the first 
Brillouin zone. It is also worth noting that this type of calcu- 



FIG. 3. (Color online) Magnetic excitations from BiFeO:! measured 
at 5 K, 140 K, 300 K, 580 K, and 680 K. The horizontal axis de- 
notes Q — Qaf- The solid lines are dispersion calculated based on 
parameters obtained from the base temperature (5 K) data. 



TABLE I. Parameters obtained from modeling the data at various 
temperatures. Note that in our spin Hamiltonian, one pair of spins 
S n and S m are counted twice, and their binding energy becomes 
H nm = t/ nm (S n ■ S m + S m ■ S n ) = 2J nm S n ■ S m , so that the 
exchange constant J used here have same definition as those used 
in R ef.flZ but should be compared to half of the J values used in 
Ref.[l3l. The units are all meV except for the ordered moment (M z ), 



which is shown as multiples of fiB ■ 
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lation is model (Hamiltonian) independent, and i?(Q) is sup- 
posed to always be 1 if the magnetic form factor used is cor- 
rect. The results from the old ionic (blue symbols) and new 
Wannier (red symbols) magnetic form factors clearly suggest 
that the new form factor provides a better description of the Q 
dependence. 
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FIG. 4. (Color online) Comparison of the original ionic Fe form 
factor and the new Wannier form factor, (a) Plotting of the ratio 

-R(Q) = j(Q^]/|f-(Qo)p ■ where ^(Q) is magnetic scattering inten- 
sity measured at wave-vector Q, while Q is the wave-vector equiv- 
alent to Q in the first Brillouin zone. Ideally this value should be 
1 (see the dashed line). The red symbols are values calculated us- 
ing the new Wannier form factor while the blue symbols using the 
old Fe ionic form factor, (b) to (d): jf(Q)| 2 plotted along different 
directions. 



Using this hybrid magnetic form factor, the ordered mo- 
ment per Fe at various temperatures is calculated and listed 
in the first row of Table [I] At T=5 K, this value M z = 4 \i B 
corresponds to Sz = 2.0, consistent with the moment ob- 
tained from other measurement :] 24 ' 25 , and close to the theoret- 
ical limit of Sz = S = 5/2, using a g— factor of 2. Frustration 
due to a positive (AFM) J2 is likely the main reason for the 
difference. We can also obtain an integrated spectral weight 
from all the spin-wave excitations, to be Sdynamic ~ 2.5 (i.e. 
M d V namic = Sdynamicg 2 ' H% = 10/x|). Combining the two 
we have a total squared moment of M 2 = AI 2 + M 2 ynamic = 
16+10 = 26 p? B , which is supposed to be S(S + l)g 2 ^i B 
to satisfy the sum rule (see Methods). This corresponds to 
S ~ 2 . 1 per Fe, in good agreement with the expected S = 5/ 2 
from Fe 3+ . On the other hand, if a simple Fe 3+ magnetic 
form factor were used, we would have obtained a Sz = 1.4 
and S = 1.4 instead, which are much lower than the expected 
values. 



Overall, our work demonstrates that the electromagnon ex- 
citations in BiFeOa can be very well modeled with a simple 
spin-wave model from a G-type AFM at low temperature. De- 
spite the cycloid modulation induced by coupling to the elec- 
tric polarization, most of the spin dynamics are barely modi- 
fied, and consistent with the behaviour of a conventional spin- 
wave. Upon heating decoherence of spin-wave excitations 
starts to occur even at just below room temperature (300 K). 
At 580 K which is well below Tjy, the coherence mode is al- 
ready almost completely destroyed. In addition, a strong hy- 
bridization between the Fe and O orbitals needs to be taken 
into account when the magnetic response is considered. Even 
as a good insulator, the electrons responsible for the magnetic 
moments are not restricted at the Fe sites, and the local mag- 
netic moment could extend well beyond the typical ionic ra- 
dius of Fe 3+ into the adjacent 2 ~ sites. 



I. METHODS 

Single crystal samples of BiFeOa are grown by floating 
zone technique at BNL. The crystal used in the measurement 
has a cylindrical shape, and weighs about 3.5 g, with a mosaic 
less than 2 degrees. Inelastic neutron scattering experiments 
are performed on the ARCS time-of-flight spectrometer at the 
Spallation Neutron Source (SNS) at Oak Ridge National lab- 
oratory. Incident energies of 60 meV and 120 meV are used 
to probe the lower and higher part of the spin excitation. 

The spin-wave calculations are based on a classical spin- 
wave model assuming a S = 5/2 G-type AFMa&. Based on 
the model, the values of J\, J r 2, and J3 can all be directly 
derived from the data if the energies of the spin-wave ex- 
citations can be accurately determined at each wave-vector. 
For instance, if we define Qi = Qaf + (0.25, 0.25, 0), 
Q2 = Qaf + (0.25,0.25,0.25),Q 3 = Qaf + (0.5,0.25,0), 
based on a simple spin-wave calculation, we have J\ = 
±(E 2 - Elf' 2 , J 2 = ±(E 2 - E 3 ), and J 3 = i(ff - 
3 Ji + 8J2). Here E\, E2, and £^3 are spin-wave energies at 
Qi. Q2. and Q 3 , respectively. 

After normalizing the data using phonon intensities which 
have been obtained from the same measurements, a sum rule 
applies: / / S(Q,haj)d 3 Qd(huj)/ f d 3 Q = S(S + 1), 

BZ BZ 

where 5(Q, faS) is the dynamic spin correlation function, and 
the integral is performed in one Brillouin zone. In our case 
/ (Q, fajj) = §<S(Q, fiu), where I(Q, ftui) is the normalized 
intensity. The factor | arises from assuming an isotropic dy- 
namic spin correlation function (i.e. S xx = S vv = S zz ), 
and averaging the polarization factor of different domains. 
The scattering intensity can be modeled with the single mode 
approximation (SMA) 27 , and the parameters that can vary 
include the exchange parameters J\, J 2, J3, ground state 
binding energies for the nearest, second nearest (face diago- 
nal), and third nearest (body diagonal) spin pairs, 2Ji(SqSi), 
2J2{SqS2), 2J 3 (SoS3), and a finite energy broadening V for 
high temperature data. 
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